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ABSTRACT 


Analytical  expressions  are  presented  which  predict  the  debonding' 
and  pull-out  lengths,  observed  in  brittle  fibre  composites.  These  charac¬ 
teristic  lengths  are  combined  with  models  of  four  toughening  mechanisms 
to  calculate  the  work  of  fracture  of  a  composite.  The  results  are  pre¬ 
sented  as  maps  showing  not  only  contours  of  toughness  but  also  the  domi¬ 
nant  toughening  micrcmechanism.  The  toughness  is  largely  determined  by 
six  material  parameters,  each  map  demonstrating  the  combined  effect  of 
changing  two  tf  these  simultaneously. 

Maps  are  presented  for  glass-fibres  in  epoxy  and  carbon-fibres  in 
epoxy.  Their  use  is  demonstrated  by  showing  the  effects  of  hygrothermal 
aging  on  the  toughness  of  the  composites. 
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Value  of  P  just  less  than  one 

Value  of  P  just  greater  than  zero 

Fibre  diameter 

Tensile  modulus  of  fibre 

General  probability  distribution  of  l 

Critical  strain  energy  release  rate  for  mode  2  fibre/matrix 
bond  failure 

Gu  divided  by  geometry  factor  a 
Shear  modulus  of  matrix 
Total  debonded  length  of  fibre 
Maximum  pull-out  length 
Mean  pull-out  length 

Average  value  of  a  distributed  length,  l 
Weibull  modulus 

Cumulative  probability  of  failure 
Cumulative  probability  of  survival 
Elastic  work  (per  fibre) 

Surface  energy  of  fibre/matrix  interface  (per  fibre) 

Pull-out  work  (per  fibre) 

Post  debond  friction  work  (per  fibre) 

General  distance  from  debond  crack 

Value  of  x  which  only  a  fraction  B  of  pull-out  lengths  can 
exceed 

Geometry  factor  from  integration  of  stress  field  around  fibre 
Constant,  but  exact  value  depends  on  author 
Surface  energy  of  the  composite 
Surface  energy  of  the  matrix 


NOMENCLATURE  (coat.) 


Failure  strain  of  fibre 

e  Failure  strain  of  matrix 

m 

n  Parameter  dependent  only  on  the  Weibull  modulus 

o  1,02  Stresses  at  which  fractions  A,  B  of  fibres  have  broken 
Debond  stress 
Average  fibre  strength 
oq  Characteristic  strength  of  fibre 

Shear  strength  of  fibre/matrix  adhesive  bond 


Frictional  shear  stress 
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1.  INTRODUCTION 

A  number  of  theories  have  been  proposed  to  account  for  the  origins 
of  toughness  of  fibrous  composites.  This  paper  describes  four  models 
based  on  the  debonding,  fracture  and  pulling-out  of  fibres.  An  equa¬ 
tion  is  derived  for  each  fracture  process  which  includes  terms  for  the 
properties  of  the  fibre,  matrix  and  interface,  together  with  the  fibre 
debond  length  and  fibre  pull-out  length.  The  fibre  debond  length  and 
pull-out  length  are  in  turn  predicted,  allowing  for  the  statistical 
distribution  of  fibre  strength.  The  paper  attempts  to  couple  our  under¬ 
standing  of  the  micromechanisms  of  fracture  with  the  models  for  each 
fracture  process,  and  thereby  predict  the  range  of  dominance  of  each 
mode  of  fracture,  and  the  corresponding  toughness.  This  is  done  by 
constructing  toughness  maps  in  which  the  axes  are,  for  example,  fibre 
strength  and  interfacial  shear  stress.  By  plotting  contours  of  tough¬ 
ness  predicted  using  the  models  for  each  micromechanism  of  fracture, 
the  effect  of  changing  the  material  parameters  on  composite  toughness 
can  be  displayed  using  the  diagrams.  As  an  example,  the  effect  of  hos¬ 
tile  environments  on  composite  fracture  behaviour  are  analysed  using 
the  maps . 

2.  THEORY  OF  FIBRE  DE BONDING 

2.1  Prediction  of  Fibre  Debondir.g  Stress 

A  number  of  workers  have  considered  the  debonding  of  fibres  in  compo¬ 
site  materials  0>2»  3,4, 5)  .  one  approach  C1*2*3*-1*)  is  .  to  consider  the  - 
stress  carried  by  the  fibre  for  which  the  corresponding  interfacial  shear 

stress  exceeds  the  shear  strength  of  the  fibre/matrix  interface,  t  .  In 

o 

this  case  the  dabonding  stress,  o^,  is  given  by: 
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The  value  of  the  dimensionless  constant  3  varies  according  to 
the  details  of  the  analysis,  but  is  independent  of  fibre  radius,  , 
in  all  cases*  Alternatively,  the  debonding  stress  may  be  determined 
by  considering  the  energetics  of  the  debonding  process.  When  the 
interfacial  shear  crack  propagates  an  increnent  6x,  the  energy 
balance  demands  that  the  surface  energy  of  the  new  interface  equals 
the  strain  energy  released  from  the  relaxed  matrix  (Fig.  1): 


2rr  (sjj.  <Sx 
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The  term  5Vm  is  the  effective  incremental  volume  of  matrix  which  is 
relaxed.  Whilst  dV^  must  be  proportional  to  the  fibre  area  and  the 
increment  6x,  its  exact  value  can  only  be  calculated  from  a  full 
stress  analysis  around  the  fibre.  For  the  purposes  of  this  work,  a 
dimensionless  constant  a  is  introduced  such  that: 


6V  =*  a  Trr^  <5x 
m 


The  debonding  stress  is  therefore  given  by: 


4E 

a  =  ( — S — IL) 

d  ar 
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(2) 


For  convenience,  let  G^/a  =»  G2,  where  G2  may  be  determined  experi¬ 
mentally. 


Equation  (1)  may  be  compared  with  the  expression  derived  by  Outwater 
and  Murphy  : 

4Ef  GII  * 
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The  dependence  on  arises  because  Outvater  and  Murphy  considered 

the  released  strain  energy  to  come  from  the  fibre. 

Note  that  the  energy-based  calculations  predict  a  debond  stress 
which  is  a  function  of  fibre  radius,  while  equation  (1)  suggests  that 
r  and  cr^  are  independent. 

Several  workers  (x*3,6)  have  ccadUcted  model  experiments  to  study 
the  debonding  of  wires  or  fibres  but  have  not  studied  the  effect  of 
fibre  size.  In  order  to  investigate  the  validity  of  equations  (1) 
and  (2)  cold 'drawn  steel  wires  of  two  different  diameters,  were  embed- 
ded  in  a  Ciba-Geigy  epoxy  system  (LY567/HY567), and  the  load  to  debond 
the  wire  measured. 

Preliminary  results  shew  the  debond  stress  is  dependent  on  fibre  dia¬ 
meter,  which  equation  (1)  fails  to  predict  (Table  1).  Equation  (2)  predicts 
a  debond  stress  ritictfor  small: large  wires  of  about  T. 9,  whereas  a  ratio  of 
2.4  was  observed  in  the  experiments.  The  difference  in  the  two  values  lies 
within  the  limits  of  experimental  error  for  the  small  number  of  speci¬ 
mens  tested.  The  energy-based  calculation  of  debonding  stress  there¬ 
fore  appears  the  more  accurate . 

In  similar  experiments,  using  bundles  of  1600  glass  fibres  of  14  ym 
diameter,  the  debonding  stress  was  estimated  to  be  0.6  (+  0.1)  GPa.  If 
the  modulus  of  the  matrix,  Em,  is  3  GPa,  then  Gz  has  a  value  of  about 
210  joules /metre2. 

2.2  Prediction  of  the  Fibre  Debond  Length 

When  the  fibre  stress  exceeds  the  debonding  stress  an  interfacial 
shear,  or  debond,  crack  propagates  along  the  fibre.  After  debonding, 
the  stress  transfer  between  fibre  and  matrix  is  no  longer  by  the  shear 
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deformation  of  the  matrix,  but  is  due  to  frictional  shear  stresses 
between  matrix  and  fibre.  This  stress  is  assumed  to  have  a  constant 
value,  xf;  in  reality  it  will  vary  due  to  the  Poisson  contraction 
of  the  fibre  within  the  matrix.  The  stress  build-up  is  therefore  assu¬ 
med  to  be  linear  with  increasing  distance  away  from  the  debond  crack 
(Figure  1). 


The  stress  in  the  fibre  is: 


4  T-  x 

°  =c°d  +  — T~  <3> 

Debonding  will  proceed  until  Q  reaches  the  tensile  strength  of  the 
fibre  end  the  total  debonded  length  is  then  given  by: 


(4) 


3 .  THEORY  OF  FIBRE  PULL-OUT 

In  the  above  discussion,  no  allowance  was  made  for  the  distribu¬ 
tion  of  flaws  in  the  fibre,  and  the  corresponding  variation  in  strength. 
The  fibre  would  therefore  be  expected  to  break  in  the  region  of 
stress;  that  is,  where  the  fibre  emerges  from  the  matrix. 


3.1  Existing  Calculation  of  Fibre  Pull-Out  Length 

Some  workers  (5»7>8)  cons£,jer  the  maximum  pull-out  length  to  be  equal 
to  half  the  critical  fibre  length.  Whilst  this  forms  an  approximate  opper 
bound  cn  it  is  not  able  to  predict  the  pull-out  length  from  materials 

properties . 
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3.2  New  Calculation  of  Fibre  Pull-Out  Length 

The  variation  of  fibre  strength  is  known  to  be  well  described  by 
the  Weibull  distribution  .  The  cumulative  probability  of  a  fibre 
surviving  a  stress  a  is  given  by: 

rr  m 

s  -  exp  (-(J-)  )  (5) 

o 

i 

and  the  cumulative  probability  of  failure  by: 

P  =  1  -  S 

A  more  refined  calculation  of  the  pull-out  length  is  made  as 
follows.  Since  the  variation  of  the  fibre  axial  stress,  given  by 
equation  (3),  is  a  linear  function  of  x,  the  pull-out  length  must 
also  be  of  a  Weibull  distribution.  Long  pull-out  lengths  occur  be¬ 
cause  there  is  a  small  probability  of  finding  a  severe  flaw  away  from  - 
in  preference  to  a  minor  flaw  near  -  the  region  of  maximum  stress. 

Figure  2  shows  how  the  variability  in  fibre  strength  leads  directly 
to  the  possibility  of  fibre  pull-out.  If  a x  is  the  stress  at  which  a 
fraction  A  1)  of  fibres  have  fractured,  and  a  2  is  such  that  a 

fraction  B  0)  have  broken,  then  a  good  estimate  of  the  maximum 

pull-out  length  possible,  is  given  by: 


The  fractions  A  and  B  may  be  regarded  as  a  measure  of  the  shortest 
and  longest  lengths  of  fibre  which  can  be  resolved  in  experiments.  This 
treatment  assumes  thac  the  whole  spectrum  of  flaws  are  repeated  in  a  length  of 
fibre  which  is  small  by  comparison  with  the  pull-out  length. 
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The  difference  in  the  stresses  and  0%  may  be  calculated 

from  equation  (3)  giving: 


(ai  -  02)  * 


4t-  l 
f  m 


The  difference  (c\  -  02)  is  a  function  of  the  variability  in  the 
strength  of  fibres.  For  a  fibre  of  uniform  strength  -  02) 

would  reduce  to  zero;  for  a  fibre  of  variable  strength,  it  may  be 
calculated  from  equation  (5)  since: 


-  -  to  S  -  - 

*  m  o 


hence: 


(°i  -  ^  *  IT  (te  B° 


The  stress  a q  is  the  characteristic  strength  of  the  fibres, 
and  is  related  to  the  mean  strength,  a^,  by: 


°  (to  2) 1/m 


Equating  (6)  and  (7)  and  substituting  for  oq  from  equation  (8) 


m  (to  2) 


a  * 

~  ,an  f  m 

TK  <B>  ‘  —J— 


Rearranging  equation  (9)  to  find  the  average  pull-out  length  &  , 
which  is  one  half  of  the  maximum  pull-out  length, 
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(10) 


where 


_  in  (A/B) 
m  (2n  2)1/m 


n  may  be  found  from  experiment;  although  A  and  B  migjit  be  esti¬ 
mated  from  the  shortest  and  longest  experimentally  resolved  pull-out 
lengths.  The  Weibull  modulus,  m,  may  be  found  from  experiment. 

Figure  3  shows  pull-out  and  debond  length  data  collected  from  90  ex¬ 
periments  on  the  fracture  of  glass-fibres  in  epoxy,  similar  to  those  des¬ 
cribed  by  Kirk,  Munro  and  Beaumont  (1C).  These  experiments  were  initially 
conducted  for  other  purposes,  but  the  results  may  be  used  to  give  an  esti¬ 
mate  of  p.  From  equations  (4)  and  (10),  the  slope  of  the  graph  is  given 


If  and  are  assumed  constant,  having  values  of  1.65 

and  0.6  GPa  respectively,  then  n  has  a  value  of  about  0.13.  The 
data-points  lying  above  the  line  are  generally  from  low  temperature  ex¬ 
periments  where  debonding  was  encouraged;  those  below  are  from  specimens 
which  had  been  exposed  to  moisture.  The  value  of  n  is  dependent  only 
on  the  reciprocal  of  the  Weibull  modulus  of  the  fibre. 

4.  WORK  OF  FRACTURE 

The  work  of  fracture  of  glass -fib  res  in  epoxy  is  thought  to  originate 
from  four  mi  crome  chan  urns  of  fracture 


4*1  Post  Debond  Friction  Energ 


The  Post  Debond  Friction  work  (12)  is  due  to  the  relative  slip  of 
the  broken  matrix  over  the  debonded  length  of  fibre,  prior  to  the 
failure  of  the  fibre.  The  work  done  per  fibre  is  given  by: 


*d  Tf  *d  fef  ‘  £m) 


The  difference  in  fibre  and  matrix  failure  strains,  (ef  -  ,  may 

be  estimated  as  the  ratio  of  fibre  strength  to  fibre  modulus  for  a 
brittle  matrix  composite, 

4.2  Elastic  Energy  in  the  Fibre 

The  fibre  stress  within  the  debonded  region  is  reduced  after 
fibre  failure#  There  is  a  corresponding  release  of  elastic  strain 
energy  from  the  fibre,  which  cannot  assist  in  crack  propagation# 

The  strain  energy,  sometimes  called  Elastic  or  Debond  energy,  was 
first  considered  by  Outwater  and  Murphy  ^  ,  and  is  given  by: 


W  -sd  f 

«  Ef  j 


I  2Tf  X| 

rf  r 


.2  K  », 


Tf  *d  Tf  af  ^d 


per  fibre 


In  order  to  simplify  the  above  expression,  it  is  often  assumed  that 
the  debond  length  is  equal  to  the  critical  length*  Equation  (11) 


then  reduces  to: 


Ttr2  ai  l . 
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i:  Since  the  debond  stress  has  been  shown  to  be  a  significant  pro- 

*  portion  of  the  fibre  strength,  the  simplification  is  invalid*  The 
full  version  of  equation  (12)  should  therefore  be  used. 

4.3  Interfacial  Surface  Energy 

The  newly  created  interface  between  fibre  and  resin  has  an  assoc¬ 
iated  surface  energy.  This  may  be  estimated  by  the  product  of  the  area 
t  of  interface  and  the  surface  energy  of  the  matrix,  y  . 


Wi  ~  ^d  Per  ^^re  (14) 

Wagner  et  al.  report  a  value  of  100  J  nT2 

4.4  The  Work  of  Fibre  Pull-Out 

The  work  to  pull  a  broken  fibre  a  distance 
socket  is  given  by: 

ird  T-  £2 

V‘> — r— 

Kelly  considered  the  average  work  of  pull-out  for  fibres 

of  length  between  zero  and  &  to  be: 

m 

l 

vU0v»  “ 

This  assumes  a  uniform  probability  of  observing  pull-out  lengthc  be¬ 
tween  the  two  limits.  Work  by  Beaumont  and  Anstice  (14)  has  shown  that 
both  pull-out  and  debond  lengths  are  well  described  by  a  Weibull  proba¬ 
bility  function. 


for  v 


&  out  of  its  matrix 
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The  calculation  of  the  average  work  of  pull-out  requires  the 
calculation  of  the  average  value  of  the  square  of  the  pull-out 
length.  Similarly,  the  average  Elastic  Energy  is  a  function  of 
V  aad  (eqn.  12). 

If  the  cumulative  probability  distribution  of  lengths  is  of 
the  form: 

S  ■  exp  <-<£-)*) 
o 


then  the  probability  of  observing  a  length  between  Z  and  (Z  +  dZ)  is 


gOO  dZ 


dS 

dZ 


dZ 


m  Z 


(m-1)  m 

exp(-(— )  )  dZ 
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The  averages  of  the  pull-out  and  debond  lengths  are  given  by  the  general 
relations: 


*  -  j  *  gOO  dJl 
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CO 

lZ  g (£)  di 
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The  required  integrals  have  been  evaluated  numerically,  and,  for 
convenience,  are  expressed  as  a  power  relationship  of  the  Weibull 
modulus.  The  expressions  shown  in  Table  2  approximate  to  the  beha¬ 
viour  of  the  average  over  the  stated  range  of  Weibull  moduli. 

The  correction  factor  for  the  pull-out  term  is  particularly  im¬ 
portant.  From  equation  (10),  the  average  length  of  pull-out  is  pro¬ 
portional  to  the  reciprocal  of  m;  the  work  of  pull-out  is  proportional 
to  the  square  of  the  pull-out  length.  The  resulting  dependence  of  W 

P 

on  m  is  therefore  given  by  the  approximation: 


/ 
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Wp  *  “72764) 
m 

In  the  case  of  the  debond  length,  the  distinction  between  7^  and 
(&)2  is  small,  and  probably  negligible  when  considering  Post  Debond 
Friction.  However,  it  may  be  important  in  the  calculation  of  the  elas¬ 
tic  work  where  the  difference  of  two  variables  is  involved. 

4.5  Other  Toughening  Mechanisms 

There  are  other  mechanisms  of  toughening  which  may  also  operate 
in  GFEP,  including  the  creation  of  fibre  and  matrix  fracture  surface 
and  their  associated  surface  energies.  These  are  small  by  comparison 
with  the  major  toughening  mechanisms. 

5.  THE  CONSTRUCTION  AND  IMPLICATIONS  OF  TOUGHNESS  MAPPING 

The  complex  interactions  of  material  parameters  cannot  easily 
be  seen  from  a  set  of  equations;  a  mapping  technique  is  therefore 
used  to  present  the  results.  A.  map  shows  the  effect  of  changing 
two  material  parameters  simultaneously ,  such  as  fibre  strength  and 
interfacial  shear  stress.  In  addition,  the  principal  toughening 
mechanism  may  be  found  for  any  combination  of  material  parameters. 

The  maps  are  constructed  by  calculating  pull-out  and  debond  _ 
lengths,  from  equations  (4)  and  (10),  for  particular  combinations 
of  material  properties.  These  results  may  then  be  used  to  calculate 
the  expected  toughness  of  the  chosen  composite,  using  equations  (11), 
(12),  (13)  and  (14),  together  with  the  results  in  Table  2. 


By  dividing  the  total  work  of  fracture  per  fibre  by  twice  the 
area  of  the  fibre  and  multiplying  by  the  fibre  volume  fraction,  the 
energy  per  unit  area  of  crack  extension,  y,  is  found. 

5.1  Exposure  of  Glass-Fibre  Reinforced  Epoxy  to  Moisture 

Contours  of  y  are  plotted  for  variation  of  fibre  strength  and 
interfacial  shear  stress  (Fig.  4),  all  other  material  parameters  are 
held  at  the  values  listed  in  Table  3.  In  addition,  to  facilitate  com¬ 
parison  with  experimental  data,  contours  of  predicted  debond  lengths 
are  plotted  for  -the  same  variables  (Fig.  5). 

These  maps  apply  to  unidirectional  glass-fibre  reinforced  epoxy 
such  as  have  been  tested  by  Anstice  and  Beaumont  (15).  As  an  example, 
they  are  used  in  this  work  to  account  for  the  fracture  behaviour  of  the 
specimens  when  exposed  to  a  100  °C/95  %  RF  environment. 

Let  us  ass  me  that  the  initial  changes  in  debond  length  are  caused 
by  resin  volume  changes  due  to  moisture  uptake,  and  further,  that  the 
fibre  strength  is  unaffected  until  cessation  of  resin  swelling  occurs 
after  3  days  exposure. 

The  trajectory  of  the  sample  on  the  x^  plane  is  shown  in 

Figures  4  and  5.  Figure  6  shows  the  predicted  and  observed  changes 
in  y.  The  prediction  follows  the  trend  well,  despite  the  crude  assump¬ 
tions  made.  The  absolute  value  of  y  is  lower  than  the  experimental 
value  because  a  small  error  in  would  cause  a  large  error  in  y, 

since  y  is  very  sensitive  to  changes  in  crfT  The  interfacial  pro¬ 
perties  may  also  be  affected,  but  there  is  ambiguity  in  the  interpreta¬ 
tion  of  changing  debond  lengths;  these  may  be  caused  by  changes  in  both 
debonding  stress  (hence  G2)  and  x^.  All  reasonable  combinations  of  G2 
and  give  approximately  equal  predicted  toughness. 
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Figure  (7)  shows  four  toughness  maps  for  glass-fibre  in  epoxy.  They 
show  ^.quantitatively*.  how  all  the  complex  interdependencies  of  material 
parameters  combine  to  affect  toughness.  They  may  be  summarised  as  follows 

i)  toughness  increases  rapidly  with  increasing  fibre  strength; 

11)  toughness  decreases  rapidly  with  increasing  frictional 

shear  stress; 

iii)  toughness  decreases  slowly  with  increasing  fibre  modulus 
and  G^ . 

In  addition,  the  toughness  is  approximately  proportional  to  the 
square  of  the  fibre  diameter,  assuming  that  fibre-strength  is  indepen¬ 
dent  of  diameter. 

5*2  Exposure  of  Carbon-Fibre  Reinforced  Epoxy  to  Moisture 

Similar  maps  have  been  produced  for  carbon— fibre  composites. 

Figures  (8,9).  These  contour  maps  are  again  compared  with  results 
from  (14).  In  the  absence  of  data  on  the  debonding  stress  of  carbon 
fibres,  a  value  of  G2  was  chosen  which  gave  good  agreement  with  the 
observed  toughness  and  pull-cut  length  for  the  first  data  point.  A 
value  of  1950  JnT2  was  used  for  G2.  This  corresponds  to  a  debond 
stress  of  2.4  GPa,  and  a  typical  debond  length  of  less  than  one  milli¬ 
metre.  Debonding  is  not  as  easily  observed  in  carbon-fibre  composites 
as  in  gxass— fibre.  This  is  partly  due  to  the  short  debond  lengths  in¬ 
volved,  and  the  different  optical  properties  of  the  fibre. 

The  changes  in  pull-out  length  and  toughness  which  were  observed 
by  Anstice  and  Beaumont  cannot  be  due  to  the  same  mechanism  as  proposed 
for  the  case  of  glass-fibre.  The  decreasing  pull-out  length  observed 
for  CFRP  exposed  to  100  °C/95°  RH  is  due  to  an  increase  in  x*.  This 
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combined  with  a  fall  in  G2  from  1950  to  ^1250  Jm  during  the 
28  days  of  the  experiment,  accounts  for  the  observed  change  in  tough¬ 
ness  (Figs.  6  and  8). 

6.  CONCLUSIONS 

ij/  The  energy  debonding  stress  of  a  fibre  appears  to  be  more  accurately 
described  by  an  energy-based  calculation. 

i£)  The  pull-out  lengths  of  a  brittle-fibre  composite  are  directly  re¬ 
lated  to  the  flaw  distribution  in  the  material. 

iii)  The  toughness  of  GFRP  and  CFRP  is  due  to  four  major  sources  which 

depend  on  six  material  constants  (o^,  E^,  d,  G2,  T£)  .  The  debond 

stress  is  also  a  slow  function  of  matrix  modulus. 

iv)  The  changes  in  the  toughness  of  composites  exposed  to  humid  environ¬ 

ments  has  been  accounted  for  using  plausible  changes  in  the  material  para¬ 
meters  <rf>  tf,  and  Gir 
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TABLE  1  —  The  dependence  of  debonding  stress 
on  wire  diameter 
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